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Abstract. We construct equilibrium networks by introducing an energy function 
depending on the degree of each node as well as the product of neighboring degrees. 
With this topological energy function, networks constitute a canonical ensemble, which 
follows the Boltzmann distribution for given temperature. It is observed that the 
system undergoes a topological phase transition from a random network to a star or 
a fully-connected network as the temperature is lowered. Both mean-field analysis 
and numerical simulations reveal strong first-order phase transitions at temperatures 
which decrease logarithmically with the system size. Quantitative discrepancies of the 
simulation results from the mean-field prediction are discussed in view of the strong 
first-order nature. 
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1. Introduction 

Networks describe interaction patterns of various complex systems. In many cases, 
interactions switch constantly and their topology changes according to external 
conditions. If we treat a network as a thermodynamic system which has an energy 
function depending on its topology and introduce temperature as the disorder strength, 
networks constitute a topological ensemble of equilibrium networks of which physical 
quantities are described by methods of statistical mechanics [H El l^- At equilibrium, 
links are reallocated to satisfy the detailed balance and ergodicity; topological phase 
transitions are expected as the temperature is varied. Since there does not exist the 
physical energy of a network associated with its topology, the functional form of the 
energy is not given a priori: It may be taken as a function of the node degree, number 
of links, or some global property of the network. If an energy function is designed for 
certain performance, one can obtain an optimized network ensemble which minimizes 
the energy via zero-temperature dynamics. On the other hand, various topologies can 
be obtained at finite temperatures, with the entropy taken into account. A number of 
networks have been constructed through the use of energy functions with a focus on 
structural transitions [H HI [5] and optimization P, [7] . In general, equihbrium networks 
constructed by a topological energy function are not complex but random at high 
temperatures; they are rather simple at low temperatures as well, making either star or 
fully-connected networks. Exceptionally, scale- free networks come out only at the critical 
temperature if a logarithmic function of the degree is introduced for the energy [1] ; this 
is logically obvious in view of that the corresponding rewiring dynamics reduces to the 
preferential attachment scheme. 

In this study, we introduce a simple function of the degree as the energy function, 
which combines two competitive terms with different strengths, and examine the 
resulting topological phase transitions. Each term favors a different ground state, so one 
may expect various network topologies depending on the temperature and the relative 
strength. 

There are five sections in this paper: Section [2] introduces the model, together with 
the energy function. In Sec. [3l the system is analyzed by means of the mean-field theory. 
Specifically, the energy and the entropy of the system are evaluated, in terms of which 
phase transitions between random, star, and fully-connected network phases are probed. 
Section H] is devoted to numerical simulations. In particular the simulation results are 
discussed, in comparison with the mean-field results. Finally, a brief summary and 
concluding remarks are given in Sec. [H 

2. Model system 

We consider a system of N nodes, some of which are connected with each other. The 
total number of links is given by M = {N/2){k), where (k) is the mean degree, i.e., 
the mean number of links per node. In this study, (k) is fixed to be 0.5 and M is kept 
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conserved during rewiring. 

Associated with the system is an energy function, which consists of two parts, with 
strengths Ji and J2, respectively. Specifically, we write the energy function in the form 



where ki is the degree of node i and {i,j) represents connected nodes. The strength of 
each part is normalized by the total number of links, so that the ground-state energy is 
extensive. This is the simplest form to describe interactions between nodes in terms of 
the degree, still affording possibility of various network topologies. 

To minimize the first part, it is advantageous for a large-degree node to be connected 
to another large-degree one. In contrast, the second part favors one node to have 
a degree as large as possible. Accordingly, the two different configurations compete 
energetically: When J2 = 0, the ground state corresponds to a fully-connected network 
for given number of links. The opposite case ( Ji = 0) yields a star network, where one 
node takes all possible links. The ratio J1/J2 thus determines the ground state to be 
either of both extremes or in-between. At high temperatures, on the other hand, the 
system is in the disordered state characterized by a random network, regardless of the 
relative strength of the two parts. Typical configurations of the three types of network 
are shown in Fig. [H 

In view of these, we expect that the system displays a topological phase transition 
between a random network and one of the two compact networks mentioned above. 
We first describe the phase transition at the mean-field level, which discloses that the 
phase transition occurs only at finite system sizes and has strong first-order nature. To 
confirm this, we carry out numerical simulations and observe discontinuity in energy, 
hysteresis, and metastability. 

3. Mean- field analysis 

We begin with the mean-field analysis of the phase transitions between the disordered 
phase (the random network) and the ordered one (either the fully-connected network or 
the star network). To describe those transitions, we introduce different order parameters, 
defined to become 0(1) in one configuration and 0(1/M) in the other. Expressed as a 
function of the order parameter and the temperature, the free energy has a maximum 
between the two extreme configurations due to the entropy contributions at finite 
temperatures. Comparing the free energies of the two configurations, we determine 
the low- and the high-temperature phases. 

3.1. Phase transition between star and fully- connected networks 

We first consider the transition between star (S) and fully-connected (F) networks and 
probe the phase boundary at which the two phases adjoin. To distinguish the two, we 




(1) 



Construction of equilibrium networks 



4 




Figure 1. Typical configurations of (a) random, (b) star, and (c) fully-connected 
networks, generated by simulations performed on a system of = 100 nodes and 
M = 25 links with (Ji, J2) = (a) (2, 1), (b) (0.5, 1), and (c) (2, 1) at temperature T = 
(a) 1.36, (b) 0.12, and (c) 0.2. 



introduce the order parameter measuring the number of nodes having the largest degree, 

= ^^^^ (2) 
n 

where it is assumed that among the n non-isolated nodes there exist only two kinds 
of nodes: Uk^^^ (= ncj)) nodes of the largest degree (fully connected to each other) and 
the remaining n{l — (p) nodes connected only to the former Uk^^^ nodes, i.e., having n(j) 
links. We consider an ensemble of a fixed total number of links, M, which are distributed 
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among n non-isolated nodes only. Then n and are related via 

nchind) — 11 



2 + n0n(l - 0) = M. (3) 
In terms of this order parameter, the energy reads approximately 

E = - ^[{n - lfnct>{n(l> - 1) + 2n{n - l){nct>f{l - 0)] 

-^[(n-l)^n0 + n(n0)2(l-0)], (4) 

which reaches the maximum between = (M + and 0^1 corresponding to a star 
and a fully-connected network, respectively. We thus compare the energy for the two 
kinds of network. In the case of a star network, we have — 1 = M, and 

Es ^ Ei<P=n-') = -^{n - 1)' - ^[{n - 1)' + (n - 1)] 

= - (Ji + J2)M - J2. (5) 
For a fully-connected network, on the other hand, we have n{n — 1) = 2M, and 

Ef = E((h=l) = ^ -n(n - if 

^ - 2JiM - 2J2V2M. (6) 

Comparison of Eqs. ([5]) and ([6]) shows that Ep becomes lower than Es for J1/J2 > 
1 — 2-^2/M + 1/M. It is thus concluded that the ground state corresponds to a 
fully-connected/star network for Ji larger/smaller than J2 in the thermodynamic limit 
(M ^ 00). 

We next consider the phase boundary at finite temperatures. To obtain the free 
energy at finite temperatures, we should take into account the entropy, defined to be 
S'(0) = Infi^. As a non-isolated node belongs to either a group of n0 nodes of the 
largest degree or the other group of remaining nodes, we write the number of distinct 
accessible configurations in the form 

The first-order transition between star and fully-connected networks is manifested by 
the barrier in the free energy -F(0) = E — TS between two competing minima. As 
addressed, the transition temperature may be determined by comparing the free energy 
for the two phases. For a star network, we have n0 = 1 and n = M + 1, thus 

Ss = \^(^Jl^^+ ln(M + 1) ^ g{a)M, (8) 

where a = 2/(A;) and g{a) = ln[a"(a — 1)^^"]. For a fully-connected network, = 1 
and the number of non-isolated nodes is given by n ~ -\/2M, which leads to the entropy 
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The free energy for each type of network thus reads 
Fs = - {Ji + J2)M -J2- Tg{a)M 

Fp = -2JiM -2J2V2M -T^^ilnM + 2\na+ 1 -ln2) 

The condition Fs = Fp at T = Tc then leads to the transition temperature 

1 



J1-J2 




(10) 



:iii 



9[a) 

below which the system makes a fully-connected network. Note that the transition exists 
only when Ji > J2 in the thermodynamic limit. This analysis does not include the true 
high-temperature phase, which is described by the random network; this is considered 
in the following section. 

3.2. Phase transition between random and fully- connected networks 

To examine the phase transition between a random (R) network and a fully-connected 
one, we define the order parameter ijj to be the average connectivity k of non-isolated 
nodes relative to the number n of such nodes: 

^ = -, (12) 
n 

where the number of total links M = nk /2 is constant, thus leading to = 2M/n^. The 
order parameter defined above conveniently characterizes the two phases, taking the 
values 1/M and 1 — (2M)~^/^ for the random and fully-connected networks, respectively. 
At the mean-field level, the energy reads 

E ^ -Jik'^ -^nk"^ = -2JiM^ -2J2^2M^. (13) 

To evaluate the entropic contribution to the free energy at finite temperatures, we write 
the number of accessible configurations 

with n = a/2M7^. With the help of Stirling's formula, the entropy is expressed as a 
function of ip: 

5 = Infi^ ^ n + nln ^ +M + Mln - 1^ , (15) 

the first two terms of which are negligible in case ip ~ 1/M. 

Note that both the energy and the entropy are monotonically decreasing functions 
of the order parameter ip. As displayed in Fig. [2], the free energy reaches a maximum 
between 1/M and 1 — (2M)~^/^ at finite temperatures. This indicates that the phase 
transition is again first-order. Taking the leading terms in the entropy, we compute the 
free energy in the two phases: 



Ff = F(^=1-(2M)~^/^) ^ -2 JiM - 2V2MJ2 - TV2M{\n V2M + 1) 
Fr = F{il)=M-^) ^ -2Ji - 2V2J2 - TM{\nM + 1). (16) 
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Figure 2. Free energy F{ijj) versus the order parameter ip, obtained from tlie mean- 
field analysis of the system of iV = 1000 nodes and M — 250 links with Ji — 2 and 
J2 = 1 at temperature T = 0.2,0.6, and 1. At temperature T higher/lower than the 
transition temperature Tc = 0.6, the free energy of a random network is lower/higher 
than that of a fully-connected network. At T — the two phases coexist and a 
first-order phase transition takes place between them. 



At the transition temperature T^, the two phases have the same free energy (see Fig. [2]). 
Equating the above two expressions, i.e., Fp{Tc) = FjiiTc), we obtain the transition 
temperature 

2J^+V27MJ2 

'~ lnM+1 ' ^ ^ 

which depends on Ji but very weakly on J2. In particular, the transition temperature 
decreases logarithmically with M, which indicates that in the thermodynamic limit 
random networks are prevalent at finite temperatures. 



3.3. Phase transition between star and random networks 

The phase transition between a star network and a random one is described by the order 
parameter [1] 

k 



M 



where k^ax is the largest degree of the network. The star-network phase is characterized 
by a large value of the order parameter, i.e., $ 1. In terms of this order parameter, 
the energy function is approximately given by 



Fr^-Lk - — ^2 



(19) 
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although it does not represent the random network phase well. To estimate roughly the 
transition temperature T^, we consider the number of accessible configurations 

This gives the leading expression of the entropy in terms of the order parameter $: 

S = -^M\nN (21) 

up to an additive constant, as reported in Ref. [1]. 

Comparing the resulting free energy function in the two phases, namely, at $ = 1 
and at $ ~ 0, we find that the former, corresponding to the star network, provides the 
global minimum of the free energy at low temperatures. The transition temperature 
below which the random network turns to a star network is given by 

J1 + J2 , . 

" IniV ' ^ ' 

which again decreases logarithmically with the size N. It is thus concluded that the 
phase transition occurs only in a finite system. 



4. Numerical results 



In this section we present results of Monte Carlo simulations performed for various 
values of J1/J2 and system size N = 100,200,400, and 1000. The average connectivity 
of the whole network is fixed to be (k) = 0.5, so that the number of total links is given 
by M = N{k)/2 = N/A. First, a random network is generated by connecting randomly 
selected two nodes M times, then the system is annealed from high temperatures via 
the standard Metropolis algorithm with randomly selected links rewired. We allow a 
danghng node to be deprived of its link and expect many isolated nodes to appear at 
the transition. Typical network configurations obtained thus are shown in Fig. [H 

Figure [3] exhibits the phase diagram of networks constructed by the energy function 
given by Eq. (1) for = 200 and 1000. It is observed in both numerical and mean- 
field results that phase boundaries depend on the system size. As predicted in mean- 
field analysis, the system undergoes a discontinuous topological phase transition from 
a random network to a compact one such as a star or a fully-connected network as the 
temperature is lowered. When Ji is smaller /greater than J2, the ground state is given 
by the star/fully-connected network. For J1/J2 ^ 1, in particular, the star-network 
phase emerges as an intermediate state, so that there occur double transitions as the 
temperature is lowered. This is more evident for small system size, which is consistent 
with the mean- field results. 

In the transition from a random to a fully-connected network of large size, it is 
observed that at the instant of sudden change of topology, the system tends to be 
trapped in a characteristic multi-star network which consists of two types of node: A 
few star nodes connected to all nodes and the remaining peripheral nodes connected 
only to the star nodes identically. Such a multi-star network appears during most 
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Figure 3. Topological phase diagram on the {Ji/T, J2/T) plane for size N —{a) 
200 and (b) 1000. Simulation results, obtained via cooling from a random network, 
are depicted by symbols, separating regions of random (R), fully-connected (F), and 
star (S) networks. Dashed lines are merely guides to the eye, and in particular those 
at low temperatures, i.e., boundaries between F and S regions, are speculative and 
just represent J1/J2 = 1, which is the minimum value to produce a fully-connected 
network as the ground state in numerical simulations. Also shown are results from the 
mean-field analysis, plotted by solid lines. 



of simulations for J1/J2 > 1 and = 1000, the largest system size considered here. 
However, it disappears if more MC steps are performed, especially on a system of smaller 
size; this indicates the multi-star network to be a metastable state. Although such a 
metastable state has also been considered in the mean-field analysis, our free energy 
function does not have any local minimum corresponding to the metastable state. 

Shown in Fig. [3] is discrepancy between numerical and mean-field results, which 
apparently grows with the system size. Here it should be noted that our numerical results 
have been obtained from cooling simulations. Considering that the phase transition is 
strongly discontinuous, accompanied with the hysteresis, we expect to obtain a higher 
transition temperature in heating simulations. With Ti and T2 denoting the transition 
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Figure 4. Hysteresis in tlie energy versus temperature curves in cooling and heating 
simulations for J1/J2 = 2.5 and N = 1000. 



temperatures in cooling and heating simulations, respectively, Fig. [3] displays that the 
annealed transition temperature Ti is lower than the mean-field transition temperature 
Tmf- In the specific case of the transition from a random network to a fully-connected 
one, Ti in numerical simulations may be estimated as follows: At temperature T = Ti, 
we assume that the slope of the free energy vanishes at ip = 1/M (corresponding to the 
random network phase) and write = dF/dip\^=i/M ~ —2JiM — a/2J2M + TiM^. This 
leads to the transition temperature in annealing 

2Ji + V2J2 



M 



(23) 



which decreases with the size much faster than the mean-field transition temperature. 
The discrepancy growing with the system size is thus explained. Similarly, at T = T2 in 
heating from a fully-connected network, we write = dF / dipl^^^i ^ — 2 JiM — J2V2M + 
T2^M/21n(av^M/2), which gives the transition temperature T2 increasing with the 
system size, 

4:V2MJi+U2 , , 

T2 = —. (24) 

lnM + 21na-ln2 ^ ' 

As a result, the hysteresis becomes more evident as the system size is increased in 

simulations. Therefore, when the system size is sufficiently large, Tmf is expected to 

locate between Ti and T2. 

We display in Fig. H] the behavior of the energy for J1/J2 = 2.5 and N = 1000, 

as the temperature is lowered (cooling from the random network) or raised (heating 

from the fully-connected network which is the ground state). Here the hysteresis is 

manifested. Note, however, that T2 obtained numerically is lower than Tmf, unlike 

the above conjecture. We presume that this inconsistency results from fiuctuations 

neglected in estimating T2. Indeed fiuctuations should assist the system to overcome 

the free energy barrier, thus lowering T2 and also suppressing its indefinite increase with 

the size. 
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5. Conclusion 

We have constructed equilibrium networks by means of a topological energy function 
which depends quadratically on the node degrees. The topological phase transitions 
between random, star, and fully-connected networks have been studied both analytically 
and numerically, through the use of mean-field and simulation methods. It has been 
observed that the system undergoes discontinuous transitions between the three types 
of network as the temperature or the interaction strength relative to the node term 
is varied. Here the transition temperatures in general decrease logarithmically as the 
system size grows. The quantitative discrepancy between the mean-field and simulation 
results is attributed to the marked hysteresis associated with the strong first-order nature 
of the transition. 
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